Effects of impurities in Spin Bose-Metal phase on a two-leg triangular strip 
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We study effects of nonmagnetic impurities in a Spin Bose-Metal (SBM) phase discovered in a 
two- leg triangular strip spin- 1/2 model with ring exchanges (D. N. Sheng et a/., arXiv:0902.4210 ). 
This phase is a quasi- ID descendant of a 2D spin liquid with spinon Fermi sea, and the present study 
aims at interpolating between the ID and 2D cases. Different types of defects can be treated as 
local energy perturbations, which we find are always relevant. As a result, a nonmagnetic impurity 
generically cuts the system into two decoupled parts. We calculate bond energy and local spin 
susceptibility near the defect, both of which can be measured in experiments. The Spin Bose-Metal 
has dominant correlations at characteristic incommensurate wavevectors that are revealed near the 
defect. Thus, the bond energy shows a static texture oscillating as a function of distance from the 
defect and decaying as a slow power law. The local spin susceptibility also oscillates and actually 
increases as a function of distance from the defect, similar to the effect found in the ID chain 
[S. Eggert and I. Affleck, Phys. Rev. Lett. 75, 934 (1995)]. We calculate the corresponding power 
law exponents for the textures as a function of one Luttinger parameter of the SBM theory. 
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I. INTRODUCTION 

There has been much interest in spin liquid phases and 
much progress has been made in our theoretical under- 
standing of these (see Ref. 1 for a review). However, only 
recently several experimental candidates have emerged. 
Among these, the triangular lattice based organic com- 
pound ft-(ET)2Cu2(CN)3 shows strong evidence of a gap- 
less spin liquid^* 3 -* 4 ^ One proposed theoretical state has 
a Fermi surface of fermionic spinons. This appears as 
a good variational stated for an appropriate spin model 
with ring exchanges and also as an appealing state in a 
slave particle studyi of the Hubbard model near the Mott 
transition, leading to a U(l) gauge theory description. 

The variational study is not sufficient to prove that a 
given state is realized in the system, and the 2D gauge 
theory does not give reliable information about the long- 
distance behavior. Driven by the need for a controlled 
theoretical access to such phases, Ref. [H considered the 
Heisenberg plus ring model on a two-leg triangular strip 
and found a ladder descendant of the 2D spin liquid in a 
wide regime of parameters; it also developed a Bosoniza- 
tion description of the quasi- ID state. 

The present work is motivated by 13 C NMR 
experiments ^ 10 i n i 12 in the organic spin liquid material 
that observed strong inhomogeneous line broadening at 
low temperatures. Theoretical Ref. [H studied effects of 
nonmagnetic impurities in the candidate spin liquid with 
spinon Fermi surface and calculated the local spin suscep- 
tibility using mean field approach. The susceptibility has 
an oscillating 2kp component decaying with a 1/x power 
law envelope. A more complete gauge theory treatment 
is expected to modify this power law j 13 i 14 i 15 but one can- 
not calculate the exponent quantitatively. 

For comparison, the ID Heisenberg chain can be 
loosely viewed as a ID version of the spinon Fermi sea 
stated and in this case the staggered component of the 



local susceptibility grows away from an impurity as x 1 ! 2 
in the limit of zero temperature and zero field. This 
was discovered by Eggert and Afflec k 17 ! 18 and is respon- 
sible for strong inhomogeneous line broadening observed 
in several ID spin-1/2 chain materials ] 19 i 2Q 

In this paper we calculate effects of nonmagnetic im- 
purities in the two-leg ladder descendant of the spin liq- 
uid using analytical approaches developed in Ref. |8|, in 
the hope of obtaining some interpolation between the ID 
chain and 2D spin liquid. We find strong enhancement of 
the 2kp components of the local susceptibility compared 
with the mean field. The susceptibility increases away 
from an impurity as ~ x l ^ 2 ~ g ^ A > x 1 / 4 , where g is one 
Luttinger parameter describing the phased and can take 
values g < 1. This is a slower increase than in the ID 
chain, but is still a dramatic effect. We also calculate 
bond textures around the defect. 



II. NON-MAGNETIC IMPURITIES IN THE 
SPIN BOSE-METAL ON THE LADDER 

The spin system resides on the two-leg triangular lad- 
der shown in Fig. [IJ which we can also view as a zigzag 
chain. Throughout we assume that the model is in the de- 
scribed descendant phase, which we will refer to as "Spin 
Bose-Metal" (SBM) following Ref. H Examples of non- 
magnetic defects are shown in Fig. [1] and are discussed in 
detail later. Generally speaking, even though there are 
different types of defects, we find that they eventually (at 
low energies) cut the system into finite sections with es- 
sentially open boundary conditions j 17 i 18 i 21 We can then 
perform analytical calculations in a semi-infinite system 
studying physical properties as a function of the distance 
from the boundary. In the following, we focus on induced 
textures in two measurable quantities - the bond energy 
and local spin susceptibility. The physics is that an impu- 
rity perturbation has components on all wavevectors and 



can directly "nucleate" the dominant bond energy cor- 
relations. The impurity also allows the uniform external 
magnetic field to couple to the dominant spin correla- 
tions, producing textures in the local susceptibility. 

Following the description in Ref. 0, there are three gap- 
less modes with the fixed-point Lagrangian density 
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Schematically, one route to this theory^ is via a Bosoniza- 
tion treatment of electrons at half- filling on the ladder, 
where we start with two bands, a = 1, 2, and assume that 
the umklapp gaps out only the overall charge mode P + 
while the other three modes P _, 0\ G , and 02a remain 
gapless. 

Ref. describes various observables in the SBM. For 
the magnetic susceptibility calculations, we will need the 
spin operator. The S z component under Bosonization is 



S z {x) 



(2) 



The most important wave vectors are Q = ±2kpi, ±2kp2, 
T(&fi + kp2) — ±tt/2, and tt. Each term can be ex- 
pressed as in Ref. [8|: 
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Throughout, we keep P + general, but it is understood to 
be pinned; details about the pinning value as well as the 
Klein factors r] aa can be found in Ref. 0. In the first line, 
the upper or lower sign corresponds to a = 1 or 2. We 
also introduce combinations a ± = (9i a ± 02a) /V2 and 
similarly for the conjugate fields ip a ±. In the last line, a 
and a' are independent numerical constants. 

When discussing non-magnetic defects and also in the 
bond energy texture calculations, we need n-th neighbor 
bond energy operator like 



B<- n \x) = S x -S x+n . 
The bosonized form can be obtained from Ref. 
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O :Nonmagnetic Impurity 
1 :Defect with Bond-Symmetry 
2, 3,4:Defect with Site-Symmetry 

FIG. 1: The top figure represents the original 2- leg triangular 
ladder model with ring exchanges, and the thick lines repre- 
sent the defects due to the impurities. The bottom figures rep- 
resent the corresponding defects in the equivalent ID model.— 
1) represents the defect symmetric with respect to a bond 
center, while 2), 3), and 4) represent defects symmetric about 
a site of the ID chain. In general, different impurities will 
lead to different fixed points. Impurity 1) will likely lead to 
a fixed point with decoupled semi-infinite systems and a non- 
magnetic cluster containing an even number of sites, while 
impurities 2), 3), and 4) will likely lead to a fixed point with 
decoupled semi-infinite systems and an effective spin formed 
by a cluster with an odd number of sites. 



where we keep only the most important wavevectors and 
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We do not show real factors in front of all terms. Here 

(n) 

7g are phases that depend on Q and the bond type: 



= nQ/2 , 



(11) 



a=l,2 



valid for Q ^ tt. Note also that since 4fc F2 = — 4/cfi 
mod 27r, there is only one independent term B^ F1 . 



A. Nonmagnetic defects treated as perturbations 

When a nonmagnetic defect is introduced at xq, we can 
treat it as a local perturbation in the Hamiltonian i 17 i 21 
Figure [1] shows some possible defects; the corresponding 



3 



perturbations are 
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Here B^ and B^ are given by Eq. © . We can charac- 
terize the defects by symmetry. In the ID chain picture, 
SH^ represents defects symmetric under inversion in a 
bond center, while 8H^ 2 ^^ are defects symmetric under 
inversion in a site. One can readily check that SH^ 2,3 ^ 
give equivalent expressions up to constant factors and, 
importantly, contain all Q modes in general. We see that 
although the defects can be characterized as two distinct 
symmetry types SH^ and 5H^ 2 \ the perturbations to 
the Hamiltonian have the same dynamical field content 
and differ only by constant phases. This is unlike the 
Bethe phase of the ID Heisenberg chain where a bond- 
symmetric perturbation contains a relevant contribution 
from a Q = tt bond operator while a site-symmetric per- 
turbation does not 

The scaling dimensions of the different contributions 
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In the Spin Bose-Metal phase we have g < 1, so the 2k Fa 
and 4/cfi terms are always relevant 0+1D perturbations, 
while the tt/2 term is relevant if g > 1/2. The relevant 
perturbations grow and one scenario is that they eventu- 
ally pin the fields at the origin. Physically this leads to 
breaking the chain into two decoupled semi-infinite sys- 
tems, which we can then study separately. The pinning 
conditions on the fields at the defect can be guessed by 
considering the most relevant perturbation and minimiz- 
ing the corresponding energy. We expect the B2k Fa and 
B^kn terms to be the dominant, which would 
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This is the case that we focus on. In Appendix [B] we will 
consider pinning conditions preferred by the term, 
which may be of interest in the borderline case g = 1. 

A comment is in order. On physical grounds, the sym- 
metry of the defect perturbation is important. For the 
case with no site inversion symmetry like the impurity 
1) in Fig. [TJ we can envision a possible outcome of the 
RG growth of the perturbation by considering a situa- 
tion where the defect bond is strong. The two spins will 



form a singlet, and if we integrate it out, we get two semi- 
infinite chains weakly coupled to each other, which under 
further RG will eventually flow to decoupled semi-infinite 
systems with pinned values of the fields at the boundary. 
We can envision more general situations where an even 
number of spins will form a strongly coupled cluster with 
a singlet ground state, and upon integrating this out we 
again have two weakly coupled semi-infinite systems. Be- 
low, we will consider a fixed point of a semi-infinite sys- 
tem and give physical calculations of the bond textures 
and the oscillating susceptibility near the boundary (im- 
purity). Turning to the case with impurities with site 
inversion symmetry like 2), 3), 4) in Fig. [TJ such reason- 
ing would give us a half-integer spin (formed by some 
effective strongly coupled cluster with an odd number of 
sites) weakly coupled to two semi-infinite systems. This 
would need to be analyzed further, which we briefly dis- 
cuss in Sec. Ill CI 



B. Physical calculations of oscillating susceptibility 
and bond textures in the fixed-point theory of 
semi-infinite chains 

From now on, we set the location of the defect to be 
the origin. We work with a semi-infinite system with 
specified boundary conditions at the origin and calculate 
the bond energy texture 

(B(x)) = (B2k Fa (x)) + (B ikF1 {x)) + (Bs(x)) .(22) 

a=l,2 

We also calculate the local spin susceptibility, which can 
be measured in Knight shift experiments. We will see 
that there are contributions that oscillate as a function 
of distance from the boundary: \{ x ) = X un% ( x ) J rX° sc ( x )] 
in fact, x° sc ( x ) dominates over x uni ( x ) an d can pro- 
duce strong inhomogeneous broadening of the NMR line- 
shapes. The local spin susceptibility at a lattice site 
i measured in a small uniform magnetic field h is 



x(i) 



dh 



h=0 



(23) 



where S* ot = Ylj $j 18 the total spin and /3 is the inverse 
temperature. Rewriting the spin operators in terms of 
bosonic fields introduced above, 



(24) 



X° sc (x)=(3 (S z osc (x) dyS z uni (y) 



where S* sc = J2q ^ x Sq and we are interested in Q = 
2km, 2k F2 , tt/2, and tt; while S* ni (y) = £ a=12 d -^M. 
Hence we define 



x osc = f3( e *Q*S z Q (x) / dyS z uni (y) + c.c.) 
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We consider the pinning Eq. (|2T|) driven by the rele- 
vant local terms B2k Fa ,B±k F1 \ in order to minimize these 



energies, the natural pinning values of 0i a (0) and #2 a (0) 
are 

cos[v / 2(9 aa (0)] = ±1 \/2(9 a(J (0) = integer x tt . (26) 

The pinning value of the field Q p - depends on the details 
such as the amplitudes and phases 7 in Eqs. ([8])- (fT0]) . 
As discussed in Appendix [Al the pinning of a at the 
origin implies stronger fluctuation of the dual field <p and 
consequently (e^^)) = 0. 

Bond energy texture is given by Eq. ([22]) . The 



(-^71-/2 (#)) term vanishes and the other contributions 
can be easily derived by applying the formulas in Ap- 
pendix [Aj 
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^sinh(^f)] 1 [f sinh(^f) 
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where a = 1,2; Aq are some amplitudes; and Sq are 
phases that depend on the pinned values at the origin 
and are ultimately determined by the details of the de- 
fect. At low temperature T — > 0, we have the following 
behavior as a function of the distance x from the open 
boundary (defect): 



(B 2 k Fa {*)) 



cos(2k Fa x + S 2 k Fa ) 
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cos(4k F1 x + S 4kF1 ) 



(29) 
(30) 



Thus, at low temperature the bond energy texture 
around the impurity reveals the correlations present in 
the system, and the physics can be viewed as a "nucle- 
ation" of the dominant "bond orders" near the defect. 
If we can tune the Luttinger parameter g, we see that 
there are two regimes: for 2/3 < g < 1 the 2k Fa terms 
dominate, while for g < 2/3 the 4:k F i dominates. 

Turning to the oscillating susceptibility, the x°/ c 2 term 
vanishes and only the X2k Fa anc ^ XV° contribute to the 
final result. Applying the formulas from Appendix [Al 
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where a = 1,2; Cq are some constant amplitudes; and Sq 
some phases absorbing all pinned field values and even- 
tually determined by the details of the defect. At low 
temperatures T — > 0, the oscillating susceptibilities at 
2k Fa and tt become 

X°2k Fa (x) - x*-icoB(2k Fa x + 5' 2kpa ) , (33) 
X% ac {x) - x°(-lf . (34) 



The envelope function in the first line satisfies x^~i > 
X4, which comes from the condition g < 1. Therefore, at 
low temperatures the oscillating susceptibility at 2k Fa ac- 
tually increases with the distance from the open end. On 
the other hand, the oscillating susceptibility at tt reaches 
a constant amplitude. 

To conclude the discussion of the semi-infinite system 
with the boundary conditions Eq. ([2T]) . we note that 
this fixed point is stable [e.g., the scaling dimension of 
^72(0) becomes 1/2 + l/(2g) > 1, so it is irrelevant]. 
The boundary spin operator has scaling dimension 1: 
e -&-> ^bound. ~ 9 x (6i a + 6 2 a) at the boundary. Know- 
ing the fixed-point theory of the semi-infinite chain, we 
can briefly discuss other situations with impurities i 17 i 22 
(For a recent review of impurity problems, see Ref. I24I .) 



C. Other Situations with Impurities 

1. Weakly coupled semi-infinite systems 

In this case, we imagine two semi-infinite chains cou- 
pled to each other at the origin. Since in each semi- 
infinite system the scaling dimension of the boundary 
spin operator is 1, the spin-spin coupling between the 
two systems is irrelevant and they will decouple at low 
energies. This is the reason why a non-magnetic impu- 
rity like 1) in Fig. [1] breaks the system into two halves 
at low energies and the physical calculations in Sec. Ill Bl 
apply generically. 

2. Spin-\ impurity coupled to a semi-infinite system 

In this case, the spin-1/2 impurity is coupled to the 
boundary spin operator which contains contributions 
from both "la" and "2a" channels, SH — \S- 1VOlV • 

( Abound., 1 + Abound., 2) * ^l^imp ' Abound., 1 ~K A2$i m p ■ 

Abound. ,2- (The "p— " sector does not enter in the impor- 
tant terms.) The couplings Ai and A2 are both marginal. 
If they are marginally irrelevant, the impurity spin will 
decouple. If one of the couplings is marginally relevant 
while the other is marginally irrelevant, the relevant cou- 
pling will grow and the impurity spin will be absorbed 
into the corresponding channel. Finally, if both of the 
couplings are relevant, since the two channels are not 
equivalent, one coupling will grow faster; a likely sce- 
nario is that the impurity spin will be absorbed into the 
dominating channel and eventually the two channels will 
decouple. 



3. Two semi-infinite systems coupled symmetrically to a 
spin- 1 impurity 

Now let us take two semi-infinite chains and couple 
them together through a spin-1/2 impurity symmetri- 
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cally. This case is also relevant for the site-symmetric 
non-magnetic impurities like impurity 2), 3), and 4) in 
Fig. [TJ The reason is that because of the site inversion 
symmetry, the non-magnetic impurity affects an even 
number of bonds which couple an odd number of spins; 
then we can imagine a strongly-coupled cluster with the 
odd number of spins, which will effectively behave as 
a half-integer spin weakly coupled to the left and right 
semi-infinite systems. 

The situation is more complex than in the previous 
subsection because we now have symmetry between the 
two semi-infinite systems, reminiscent of the 2-channel 
Kondo problem. We can imagine the following possibil- 
ities. When all couplings are marginally irrelevant, the 
impurity spin and the two semi-infinite systems will de- 
couple at low energies (and the physical calculations of 
textures in Sec. Ill Bl are valid in this case). Suppose now 
we have marginally relevant couplings and the dominant 
growth is for the channels la in the two semi-infinite 
systems. One is tempted to speculate about the possi- 
bility of "healing" the channels la across the impurity, 
while the channels 2a remain open. However, it is likely 
that this is not a stable fixed point in the presence of 
the allowed terms in the Hamiltonian coming from the 
microscopic ladder system. While the eventual outcome 
is not clear and depends on details, on physical grounds 
we again expect arriving at some stage at a fixed point 
with some odd number of spins forming a half-integer 
spin that is decoupled from two semi-infinite systems. 



III. CONCLUSIONS 

To summarize, following the theoretical description 8 of 
the Spin Bose-Metal phase in the triangular strip spin- 
1/2 model with ring exchanges, we discussed the effects 
due to different types of impurities. The defects can have 
additional bond or site symmetry in the ID zigzag chain 
language. We first treated the defects as local perturba- 
tions in the Hamiltonian and saw that all types produce 
relevant perturbations, eventually breaking the system 
into two halves and a separate decoupled cluster of spins. 
In the bond-symmetric case (or more general cases with 
no symmetries) the decoupled cluster is likely to be non- 
magnetic, while in the site-symmetric case it has half- 
integer spin, and the details of such fixed points depend 
on the microscopic details j 17 i 22 This analysis also moti- 
vated appropriate boundary conditions for pinning the 
fields in the fixed point theory for the semi-infinite sys- 
tems. 

For such a semi-infinite chain, we calculated the bond 
energy texture near the boundary and found power 
law decays Eqs. (|2T]) -(|28 ]) of the oscillating components 
at wavevectors 2k Fa, 4/cfi- The dominant power law 
switches from the 2k Fa to the when the Luttinger 
parameter g drops below 2/3. We suggest that charac- 
terizing such bond textures in numerical studies 8 could 
be useful for determining the Luttinger parameter g of 



the SBM theory. 

We also calculated the oscillating susceptibilities at 
2kp a and 7r, Eqs. ([3T]) - ([32|) . which behave differently at 
low temperatures. The susceptibilities at 2k Fa actually 
increase with the distance from the boundary in the limit 
of zero temperature (and zero field), while the suscepti- 
bility at 7r becomes distance-independent. The rate of in- 
crease at 2kp a is slower than in the ID chain^ but would 
still produce strong NMR line broadening at low temper- 
atures. Of course, this is the result for the long-distance 
behavior along the ID direction. If we are thinking about 
the 2D spin liquid, we would likely expect a power-law 
decay away from an impurity j 13 i 14 i 15 Nevertheless, the 
persistence of the oscillating susceptibilities on the quas- 
1D ladders suggests that in the 2D case the decay may 
be slow and also produce significant inhomogeneous line 
broadening. 
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APPENDIX A: ONE MODE THEORY ON A 
SEMI-INFINITE CHAIN 

For the simplest case / 7 i 18 i 23 consider a one mode the- 
ory on a semi-infinite chain with pinned value at the ori- 
gin, #(0, r) = pinned = 6$. The action is 



poo pf3 

S = / dx dr 

Jo Jo 2tt# 



v(d x 0) 2 + -(d r 0) 



. (Al) 



The correlation functions needed in this paper are 



iuQ(x,r) 



[f sinh(^)]V 
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7T / (3V 
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Pv [ ^ sinh( ^ )] ^ 



(A2) 

iue(x,r)y A3 } 

, (A4) 



o . 



(A5) 



Here u is a parameter depending on which quantity is 
being measured and A is some real constant. The {e %u0 ) 
is non-zero because of the pinning at x = and decays as 
a power-law away from the origin at T = 0. On the other 
hand, the conjugate field ip fluctuates more strongly than 
in the bulk and (e zuip ) = everywhere. 

We can similarly consider a one mode theory with the 
dual field pinned at the origin, <^(0,r) = pinned = cpo. 
It is convenient to work with the action 



dx dr — 
Jo ^ 



v 



(A6) 
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The correlation function needed is 

lUip(x,T) \ 



[f sinh(^)]%T 



(A7) 



where A is some constant. 



APPENDIX B: CALCULATIONS IN A FIXED 
POINT THEORY OF A SEMI-INFINITE SYSTEM 
WITH PINNED 6W(0), ^-(0), AND <p p -(0) 

Here we consider the theory Eq. (pp) on a semi-infinite 
chain with boundary conditions 



Pin^+(0), (Ax-(O), (p p -(0) . 



(Bl) 



This can arise if we minimize the B n /2(0) perturbation 
instead of the B2k Fa (0) and £?4*; F1 (0), see Eqs. (|8lfT0|). 
Coming from the microscopic ladder spin system in the 
SBM phase with g < 1, this fixed point is unstable to 
the allowed i?2fc Fa (0) terms. Nevertheless, it can be of 
interest in the special case with g = 1, which is realized, 
e.g., by the Gutzwiller wavefunctions or at phase transi- 
tions out of the SBM. 8 The calculations of the physical 
textures are simple and we summarize these below. 

Because of the pinning Eq. (|Bip , the dual fields (pa+, 
6 a -, and p - fluctuate more strongly. The only non- 
vanishing term in the bond energy texture is 



(B(x)} = (B n/2 (x)) 



(B2) 



If v\ = V2, then "cr+" and "cr— " variables decouple and 
we can apply the formulas in Appendix [Aj In the general 



case vi V2, the calculations are more demanding but 
the result is simple: 



A*/ 2 COS (fx + 6 w/2 ) 



iv£ s inh(^|) ^sinh(^) 

7T v V r fj 'IT v Vfj ' 



,(B3) 



where ^- = ^ (t?T ^ 1^) ' Ar/2 is some amplitude and 
5^/2 is a constant phase. The pinned values at the 
origin as well as the Klein numbers enter in the same 
way as they enter the assumed minimization of 1^/2(0), 
Eq. ([TO]) , so the final result depends only on the physical 
details of this term at the origin. In the limit T — > 0, 



(B«/2(x)) 



cos (fx + 5^/2) 



(B4) 



As for the oscillating susceptibility, similarly to the 
bond energy texture, only the Q = tt/2 term contributes 
to the final result. Again, for v\ = V2 we can apply the 
formulas in Appendix [Aj while in the general case we get 



Xtt/2 



C • x • cos(|x + S' n , 2 ) 



sinh(^|) M sin h(2zifi) 



(B5) 



where is some amplitude and 5' n , 2 is a constant 

phase absorbing all pinned values and eventually deter- 
mined by the details of the defining energy £^72(0). In 
the low temperature limit T — > 0, 



X%%~x* 4 *cos(^ + £; /2 ) . 



(B6) 
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